Vortex Lattice Inhomogeneity in Spatially Inhomogeneous Superfluids 
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A trapped degenerate Bose gas exhibits superfluidity with spatially nonuniform superfluid den- 
sity. We show that the vortex distribution in such a highly inhomogeneous rotating superfluid is 
nevertheless nearly uniform. The inhomogeneity in vortex density, which diminishes in the rapid- 
rotation limit, is driven by the discrete way vortices impart angular momentum to the superfluid. 
This effect favors highest vortex density in regions where the superfluid density is most uniform 
(e.g., the center of a harmonically trapped gas). A striking consequence of this is that the boson 
velocity deviates from a rigid-body form exhibiting a radial-shear flow past the vortex lattice. 



It has long been understood that a superfluid can only 
rotate by nucleating quantized vortices P, Q , which in 
turn control most of its macroscopic properties. In recent 
years, rapid progress in the field of confined degenerate 
Bose gases has led to the experimental realization of largc 
vortex lattices H II II H 

One of the most striking features of these arrays is thcir 
apparent uniformity, despite the strong spatial variation 
of the local superfluid density imposed by the trap. These 
observations cannot be simply explained by the strong 
vortex interaction, nor by appealing to the imposed rigid- 
body rotation. Based on purely energetic considerations, 
one would expect a highly nonuniform vortex distribution 
that is suppressed at the center of the trap where the su- 
perfluid density (and therefore kinetic energy cost), as 
well as vortex repulsion, are largestQ- Despite some at- 
tempts to understand this uniformity @, which has also 
been observed in simulations [lOj, no clear physical ex- 
planation has appeared in the literature. 

In this paper we present a theory of vortices in a con- 
fined, spatially inhomogeneous rotating superfluid. We 
provide a simple physical explanation for, and compute 
corrections to, a uniform vortex array. Our main result 
is the vortex density n v (r) which, in the simplest case of 
a harmonie trap with a strong uniaxial anisotropy and 
within the Thomas- Fermi (TF) limit, is given by 
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with R the TF radius, u> = Om/U the rescaled rota- 
tional velocity f2, m the boson mass, and £ the coherence 
length [ïï| . 

The nearly uniform vortex distribution n v ra u/n is 
a consequence of a balance between spatial variations of 
the kinetic energy per vortex and the vortex chemical 
potential, both of which scale with the local superfluid 
density, p s (r). While it is energetically costlier to posi- 
tion vortices in a region where p s (r) is high (the center 
of the trap), the vortex chemical potential (controlled by 
p s (r)u>) is also high there, compensating and leading to 
an approximately uniform vortex density. 

A spatially-dependent correction to n v (r) in Eq. |JT|| 
arises from vortex discreteness and the related inability 
of the vortex state to locally reproduce uniform vorticity 



corresponding to rigid-body rotation. In an inhomoge- 
neous condensate, the associated kinetic energy-density 
cost is spatially dependent, and is lowered by a nonuni- 
form vortex distribution. The reduction in n v {r) scales 
with V 2 lnp s , i.e., it is smallest where the condensate is 
most uniform, and leads to the strongest vortex density 
suppression away from the center of the trap. The correc- 
tion vanishes in the uniform condensate (R — ► oo) and, 
(within the London approximation) dense vortex (fast 
rotation, w£ ~ 1) limits, in which condensate inhomo- 
geneity and vortex discretness are (seemingly) unimpor- 
tant 

An immediate interesting consequence of the radial 
vortex lattice distortion is that the corresponding az- 
imuthal superfluid velocity v s (r) deviates from a rigid- 
body form, exhibiting radial-shear flow. We expect that 
this will induce an azimuthal shear distortion of the lat- 
tice, with chirality set by the sense of the imposed rota- 
tion Below we sketch the derivation of these results. 

A rotating superfluid is most easily analyzed in the 
frame in which the boundary conditions (i.e. the prover- 
bial bucket) are stationary. For experiments on trapped 
Bose gases, this is the frame rotating with frequency f2 
(in which the normal fluid is stationary 0). For sim- 
plicity we focus on a trap with a high degree of uniaxial 
anisotropy, which reduces the problem to two-dimensions 
perpendicular to the (z-) axis of rotation. Deep in the 
superfluid state the London description, which focuses on 
the superfluid phase 9 degree of freedom, is sufficiënt and 
is represented by the energy 



E= — / d 2 rp s (r)[(Vef - 2u(z x r) • V0], 
2m J 
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where within the TF approximation p s (r) ra (p—V(r))/g, 
with /i the boson chemical potential, g the s-wave scatter- 
ing potential and V(r) = Vr{r) — ^mfl 2 r 2 a combination 
of the trapping (Vr{r)) and "centrifugal" potentials. 

Before addressing the many vortex problem, we present 
the solution for a single vortex in an inhomo gen eous su- 
perfluid, which, despite a number of studies |l4l ITEl ITfiL 
IÏtI IHjj . has eluded a complete solution. For a vortex 
located at r off the trap center, the superfluid veloc- 
ity v s = (?i/m)VÖ is determined by the Euler-Lagrange 
(EL) equation for Eq. @ (expressing boson number con- 
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servation) V- (p s (r) Vö) = 0, subject to the vorticity con- 
straint V x Vö = 27r<5^ 2 ^(r — r ), and a boundary condi- 
tion of a vanishing superflow transverse to the boundary. 
In contrast to uniform systems found in condcnsed mat- 
ter experiments (c.g., helium in a bucket), in an atomic 
trap p s (r) vanishes at the boundary thereby automati- 
cally satisfying the boundary condition. The solution can 
then be expressed as 9 = 9 V + 9 a with V9 V = z * r < ^~*y 
the usual (uniform-superfluid) vortex form ensuring the 
topological vorticity constraint, and the single- valued an- 
alytic phase 9 a satisfying the EL equation: 



V Ps ■ V9 a + Ps V 2 9 a = -V Ps - V6 



(3) 



We reserve the full analysis of Eq. © to Ref . 0, focusing 
here on main results. Near ro, we find, in agrccmcnt with 
Refs.G3H 
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Far away from the vortex, v° = (h/m)'V9 a (r) w 0, van- 
ishing like a dipole field with negative and positive vor- 
tices at ro and at the center of the trap (r = 0), respec- 
tively. Because for a typical trap Vp s (ro) cc — ?o, and the 
logarithm is negative for r — > ro, the analytic distortion 
in the superfluid velocity is perpendicular to ro and leads 
to a superflow that is no longer purely azimuthal around 
an off-axis vortex. In agreement with local mass conser- 
vation, v s is smaller on the trap-center side of the vortex 
(where p s (r) is larger) and larger on its outside (where 
p s (r) is smaller). A refinement of experiments 0, E3 
that have demonstrated the ability to measure the phase 
variation 9(4>) around a vortex should allow a direct de- 
tection of the superflow distortions predicted here. 

In an ideal superfluid and in the absence of other forces, 
a vortex moves with the local superfluid velocity. Therc- 
fore, the above result has interesting implications for vor- 
tex dynamics. Namely, since v" (r) is finite at the center 
of the vortex (where |r — r 1 w £), we find a remarkable 
result: a single vortex at radius r, without (sustained) ex- 
ternally imposed rotation will precess about the trap cen- 
ter at a frequency w « (h/2mr)(d r p s (r)/2p s (r)) ln « 
(^/2mi? 2 )ln|, that away from the condensate edges is 
roughly independent of r. Such vortex precession has in 
fact been seen experimentally |2l| . with a quality factor 
of order 10. The implication of the superflow distortion 
v"(r) is even richer for dynamics of a pair of vortices. We 
predict that two same-charge vortices will orbit their cen- 
ter of charge, that will in turn precess about the center of 
the trap, with two frequencies determined by vortex sep - 
aration and location relative to the axis of the trap [ïjj . 

We now turn to the many-vortex problem, with the 
goal of computing the vortex spatial distribution in an 
inhomogeneous rotating superfluid. The total v s , mea- 
sured in the laboratory frame, due to an array of N vor- 
tices is the sum of the contributions from each vortex: 
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with vortex positions r^ static in the frame of the nor- 
mal component. In the above, we have neglected v°(r), 
as it has a subdominant effect in the many-vortex prob- 
lem |T^ . The corresponding vortex density is given by 
n v (r) = (2^)-^ x V9 = £f 5 2 {v - n). 

For large f2 the vortex state is dense, and we can ne- 
glect the discrete vortex nature [embodied by Eq. (JSJ] and 
approximate v s (r) and n v (r) by arbitrary smooth func- 
tions. Expressing E in Eq. (J2J) in terms of v s (r) (and 
dropping a constant), we have 
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which is clcarly minimized by the rigid-body solution 
v s = Qz x r corresponding to a uniform vortex density 
n V Q = uj/t: = mQ/nh. 

Away from this classical rapid-rotation limit, vortex 
discretness begins to matter and the above solution 
clearly breaks down, as v s (r) diverges as 1 / j r — r ^ j near 
each vortex at r^; it thus strongly deviates from rigid- 
body flow. In this regime, where a superfluid exhibits its 
locally irrotational quantum nature, the summation in 
Eq. (JSJ can no longer be replaced by an integration, and 
the minimization of E must be done directly over the r.; , 
rather than over a field v s (r). 

For a uniform infinite condensate the problem was 
solved long ago by Tkachenko j^, who found that the 
solution is a hexagonal lattice characterized by the vortex 
density ü v q. Even in the uniform case, for a finite system 
vortex discreteness manifests itself in the lower-critical 
rotational velocity Q c i (H/mR 2 ) ln j- below which no 
rotation is supported by the condensate. 

To analyze an inhomogeneous condensate, it is essen- 
tial to faithfully incorporate vortex discreteness in treat- 
ing the sum in Eq. © . For r near a vortex located at Tj , 
the flow is dominated by a diverging contribution from 
the j-th vortex, with other vortices giving a subleading 
correction that is smooth for r — * r j . With this observa- 
tion v s (rj + Sr) near r^ is well-approximated by 



v s (rj + Sr) 



h z x Sr 
m Sr 2 



v s(rj), 



with the smooth superflow v s (r) 



v s (rj) ~ - / rfVn„(r') ; , 
m J (rj — r'y 



(7) 



(8) 



due to all other vortices, expressed through a coarse- 
grained vortex density n v (r) — n v o(l— V-u(r)), or equiv- 
alently vortex displacement u(r) that are to be deter- 
mined fïjl . Interestingly, and not unlike the special treat- 
ment of the discrete (k = 0) ground-state BEC mode, 
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to retain vortex discreteness we extracted the dominant 
contribution to the flow around the j-th vortex and then 
safely replaced the rest of the sum by an integral over the 
coordinates of the other vortices. By taking a curl with 
respect to Tj and solving for v s (r), Eq. JÏÏJ) gives 



v s (r) = Q[z x r — 2z x u(r)], 



(9) 



describing the deviation of the mean-field flow from the 
rigid-body velocity due to the distortion u(r) from a uni- 
form vortex lattice. 

To compute the optimum vortex density n v (r), we ex- 
press the total energy E of a vortex array [i.e., Eq. J5J] 
as a sum over lattice cells[23|, with each associated 
with a single vortex. Using Eqs. (|7I9|I for v s inside a 
cell, making a circular-cell (of area l/n„(r :) )) approxi- 
mation, and assuming p s (r) does not vary appreciably 
on the scale of the vortex spacing, we compute the en- 
ergy per cell. The remaining sum over cells can be eas- 
ily done, inoccuously approximating it by an integral, 
i- e -j J2i ~* ƒ d 2 rn v (r) ss ƒ d 2 rn v (r). Thus we obtain 



E~ / d 2 rp e (r) w(l - V-u) ln 
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£ 2 w(l- V-u) 



Alü 2 u 2 
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where we have discarded u-independent terms. 

Minimizing £7[u(r)] over u(r), and re-expressing the 
solution in terms of n v (r) , we finally find 



nv(r) = - + —V 

tt 8ir \ps(r) 



' V[p s (r) ln 



7r£ 2 n„(r) 
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that gives a local vortex density in a rotated trapped 
superfluid, characterized by a rotation rate f2 and local 
superfluid density p s (r), i n ag reement with recent lowest 
Landau level results l25l l26l l27| . For a smoothly- 
varying p s (r), we expect that vortices will have a locally 
hexagonal lattice structure [22^. but with a lattice pa- 
rameter a(r) = (2/V3n„(r)) 1 / 2 that varies with radius. 
Equation consists of a uniform contribution corre- 
sponding to rigid-body rotation of the superfluid and a 
correction that depends crucially on the p s (r) profile and 
therefore on the shape of the trapping potential. As is 
clear from its derivation, this spatial variation of the vor- 
tex density arises from an interplay of the nonuniform 
trapping potential and vortex discreteness, the latter a 
fundamentally quantum-mechanical effect. 

For a rotating condensate of size R, trapped in a 
smooth concave potential, p s (r) varies on the scale R, 
leading to a negative correction (i.e. the second term in 
Eq. I|11|0 to the rigid-body vortex distribution, that van- 
ishes in the thermodynamic limit as l/R 2 . Therefore, 
(as depicted in Fig^l we predict, generically, a vortex 
distribution that is denser at the center of the trap and 
falls off towards the condensate edge as V (Vp s (r) / p s (r)) . 
Within our London approximation, this nonuniformity 
also vanishes in the fast rotation (w£ 2 ~ 1), dense- vortex 
limit, in which vortex cores overlap and their discreteness 
is unimportant |l2j . 
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FIG. 1: (Color online) Main: Vortex density as a function 
of radius (normalized to n„o, dashed line) for these Q and Ti 
values (labelled by f2): f2 = .86f2T and R = 49um; f2 = .57f2r 
and R = 31/im; points are experimental values for the Q = 
.86f2T case adapted from Ref. l28l Inset: Hexagonal lattice 
parameter a/R as a function of radius for the f2 = .86f2T 
case; points correspond to same data as in main. 



Recalling that a uniform vortex distribution corre- 
sponds to a rigid-body rotational superflow, v s (r) = 
f21 x r, the vortex lattice distortion predicted in Eq. (TT} 
has an immediate remarkable consequence, namely a ra- 
dial shear of the superfluid velocity, with the average ro- 
tational frequency descreasing with radius. As a result, 
unlikc a uniform lattice, a radially distorted vortex lat- 
tice rotating at rate f2 cannot be stationary with respect 
to the surrounding shearing superfluid. Symmetry argu- 
ments then suggest that for a nonideal superfluid, a vor- 
tex lattice should exhibit a chiral azimuthal radial-shear 
distortion in the direction opposite to the fluid flow |13| . 

We can use Eq. (|ll|l to predict the radial vortex density 
for a rotating harmonically trapped BEC with trap fre- 
quency f2*r. Deep below T c for a harmonie trap and away 
from the condensate edges we expect p s (r) to be well ap- 
proximated by the TF expression p s (r) ~ p®(l — r 2 / R 2 ). 
Here, R~ 2 = ■m(Q 2 ^—fl 2 )/2p is reduced by the applied ro- 
tation; however, we shall regard it as an experimentally 
given parameter. Together with a mild approximation 
n v (r) sa lu/tt in the argument of the logarithm, we imme- 
diately obtain Eq. (|TJ, which in Fig. [I]we plot for 87 Rb 
using realistic 28] parameters of f2 = .86f2T, R = AQpm 
(top curve), and f2 = ,57f2T ; R = 3lpm (bottom curve) 
along with data from Ref. 1281 for the farmer case. The in- 
set shows a for the f2 = .86f2r case. Here, £ = y/h/mQ,T 
is the TF value and f2-r = 52s _1 . This agreement 
with experimental data is achieved with no adjustable 
parameters. 

A combination of magnetic and optical trapping allows 
an unparalleled degree of control over the single-body 
potential seen by the atoms. A wide range of experi- 
mentally accessible trapping potentials allows stringent 
tests of our predictions. Consider, for example, a som- 
brero potential, V(r), consisting of a shallow overall trap 
with frequency f2 j 1 and a repulsive Gaussian of width £, 
with £ <C R- The superfluid density will display a cor- 
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responding "dip" near the trap center, which within the 
TF approximation is given by 

Ps{r) = Pso - Psi cxp(-r 2 /2^ 2 ), (12) 

with p S Q (> p a \) arising from the shallow part of the 
trap and therefore approximately constant on the scale t. 
Using this p s (r) inside Eq. we find a vortex density 

n v /n v o 
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FIG. 2: (Color online) Main: Fhe solid line depicts the vor- 
tex density (normalized to n„o; dashed line) as a function of 
radius for the condensate profile Eq. 1121 . Inset: Locations 
of vortices (with lines to guide the eye). 
profile n v (r), that is, interestingly, non-monotonic and 
near the center of the trap exceeds the asymptotic rigid- 
body value[l9|. We plot n v (r) for a large, slowly rotating 
condensate (R = 500/im, SI — .2s _1 , p s i/p s o = 0.67, and 



l = .3R) in Fig. The growth in n v (r) near the cen- 
ter of the trap manifests itself as contracted inner vortex 
rings, as shown in the inset. Unfortunately, current ex- 
periments, which work with relatively small condensates 
and vortex lattices, quickly run out of length scales to 
quantitatively test details of a spatially varying vortex 
distribution. We hope, however, that this qualitatively 
distinctive vortex response to a tunable trap potential 
will be observable in the next generation experiments. 

To conclude, we have studied the vortex spatial distri- 
bution in a rotating trapped superfluid and showed that, 
consistent with experiments and despite a strongly inho- 
mogeneous condensate, in the limit of high rotation rates 
and a large condensate, vortex density is nevertheless 
nearly uniform. We have computed the leading spatially- 
dependent correction to n„(r), which arises from an in- 
terplay of an inhomogeneous trap potential and vortex 
discreteness, and showed that generically it leads to a 
vortex density that is largest at the center of the trap, 
in striking contrast to simple energetic expectations. We 
hope that our predictions, which are quantitatively and 
qualitatively consistent with recent experiments |28|, will 
stimulate more detailed experimental work on inhomoge- 
neous vortex states. 
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dington, E. Cornell, P. Engels, A. Fetter and V. 
Schweikhard, as well as support from NSF DMR-0321848 
and the Packard Foundation. 
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